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The dynamical growth rate of bubbles nucleating in rclativistic plasma in thermal 
first-order phase transitions is analyzed. The framework is a hydrodynamical model 
which consists of relativistic fluid and an order parameter field. The results of 
analytical approximations and numerical simulations coincide well. 



a To appear in the Proceedings of Strong and Electroweak Matter '97, Eger, Hungary, 
21-25 May 1997. 



In thermal systems first-order phase transitions normally proceed via nu- 
cleation of bubbles of the new phase. Bubbles larger than a certain critical 
size begin to grow, whereas smaller bubbles will shrink. Langer's formulaBH 
for the nucleation rate of bubbles of the new phase is given by 

r = e W (1) 

where T is the probability of nucleation per volume and time, k a dynamical 
and VLq a statistical prefactor, AF the free energy difference of the system with 
and without the nucleating bubble, and T the temperature. The purpose of 
this study is to investigate the dynamical prefactor k. 

Let the phase transition be driven by an order parameter field <p(x). Initial 
growth rate of perturbations around an extremum configuration is given by the 
coefficient k: 

S(f> oc e Kt . (2) 

Similar relations hold for perturbations of energy density and fluid velocity. A 
parallel definition for the growth coefficient n can be expressed in terms of the 
radius R(t) of an expanding spherical bubble of the new phase, 

^««[i?(t)-i? cr ], (3) 

where i? cr is the radius of the critical or extremum bubble. This equation is 
valid when R(t) i=s R CT . In order for the bubble to grow, the initial radius R(0) 
must be slightly larger than that of a critical bubble. 

In literature there are at least two calculations of the growth rate k for 
relativistic plasma. Csernai and Kapusta write the free energy density as a 
functional of their order parameter which is just the usual internal energy 
densityd Ruggeri and Friedman make the approximation that the interface 
is infinitely thin and do not need any order parameters These two methods 
produce results for k which disagree even qualitatively with each other. 

Let us now introduce a hydrodynamical model which enables both an 
analytical and a numerical determination of the initial growth rate k. The 
modelu consists of an order parameter field (f>, and perfect fluid which describes 
the other degrees of freedom. There are three basic locally varying quantities, 
namely </>(x), fluid four- velocity ^(x), and temperature T(x). Due to the 
small value of the baryon asymmetry all the conserved particle numbers can 
for the present purposes well be approximated to be zero in the early Universe. 
Furthermore, the expansion of the Universe can be neglected, since the whole 
period of nucleation, yet alone the initial growth of bubbles, is an extremely 
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rapid process. Thus the equations of motion can be written in the form 
f d T^ v = 

( 0V + £y(0,T) = V> (4) 

where T^ v = T^ v {u a (x), dg<j)(x), <j)(x), T(x)} is the energy-momentum tensor, 
V the potential energy density for the order parameter field, T the tempera- 
ture, and t] a phenomenological friction parameter (not to be confused with 
shear viscosity). The upper equation is the conservation law of total energy- 
momentum. The lower equation tells how energy is transported between the 
order parameter and the fluid through the dissipative term, proportional to rj. 
The same term is also responsible for the creation of entropy. In the limit of 
vanishing fluid velocities the lower equation gives the simple dissipative equa- 
tion 

d0_ I 6S3W (5) 
dt 77 6cj> 

where S3 [<j>\ is the usual high-temperature three-dimensional action (equalling 
free energy at extrema). 

It is clear that a hydrodynamical description of the system cannot be 
complete. It is only valid at scales which are longer than particle mean free 
paths or interaction times. The solutions to Eqs. (0) are smoothly behaving 
fields, whereas in reality fields have strong thermal fluctuations on short scales. 
To incorporate thermal fluctuations a Langevin-type equation with a noise 
term would be needed. Within a purely hydrodynamical model questions for 
example about damping effects in plasma cannot be answered. 

The model in Eqs. (^) can be applied to describe both electroweak and 
QCD phase transition. In electroweak theory the order parameter is obviously 
identified as the Higgs field. The true coupling term is more complicated than 
the frictional ry-term of this model. But that effect should not be significant, 
as long as the detailed internal structure of the interface is not being discussed 
(something which a hydrodynamical model cannot accurately determine any- 
way). The value of the friction parameter rj can be fixed by comparing with 
microscopic calculationsQ In the case of QCD the order parameter cannot be 
identified with a physical particle. However, one can still employ the model as 
a purely phenomenological construction. Interaction length or time of QCD 
sets the scale for the friction coefficient r\. A further complication is the macro- 
scopic mean free path of neutrinos, but luckily the hydrodynamic energy flux is 
in normal cases clearly superior compared with that carried by the neutrinos. 

For solving Eqs. ([|) the potential V(<p, T) must be known. Here the usual 
quartic potential has been employed. By fixing the parameters of it in a suit- 
able manner, the desired values for the surface tension and latent heat of the 
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transition will be reproduced. The coordinate system is spherically symmetric 
1+3 dimensional space-timeQ 

In order to create the initial configuration the critical bubble solution must 
be known with high accuracy. Going closer to the thin- wall limit, that is, using 
larger critical bubbles, has the advantage that numerical errors in the deter- 
mination of growth coefficient k decrease. But in this limit the field equation 
cannot be directly integrated numerically to produce the critical bubble. In- 
stead, the following Ansatz is used: 

W'I-t{ 1 -'«*(^)}' (6) 

Here i? cr , £ cr are unknown parameters, and m i n is the position of the new 
minimum of the potential. In this two-dimensional subspace set by the Ansatz 
the extremum of the action becomes saddle point of an ordinary function, 
Si[(j)A\ = S3a(Rci, £cr)- This saddle point can then be located numerically. 
For larger bubbles this method produces quite accurate results, which was 
actually unexpected. Fixing R cr naively by Laplace's relation, i? cr = 2a/ Ap, 
where a is surface tension and Ap pressure difference between the two phases, 
would lead to huge inaccuracies. However, the correct value of_i? cr can be 
found directly by employing curvature-dependent surface tension!!, a(R). 

Analytically the initial growth rate k can be determined as follows. Expand 
the low- velocity dissipation equation (||) around the critical bubble </> by making 
the substitution <f)(t, x) = 4>(r)+Lp(t 7 x). The resulting equation for fluctuations 
is 

d^_ is 2 sS] 



dt rj S 2 4> 



V (7) 



Next insert the unstable growth mode, Eq. (g). The thin- wall result for the 
negative eigenvalue of the fluctuation operator, A = 2/R 2 r , has been observed 
to hold surprisingly well in the general case, tooEl This gives the approximate 
solution for the initial growth rate in the hydrodynamical model: 

o 

(8) 



Radius of the critical bubble, R cr , depends on one hand on the cooling rate 
of the system — in cosmology on the strength of gravitational interaction — and 
on the other hand on the thermodynamical properties of the phase transition, 
especially on the values of latent heat and surface tension. 

A more straightforward method is to integrate the dissipation equation (^|) 
directly in the thin-wall limitlI3 The result for the initial acceleration of the 
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Figure 1: Numerical determination of the initial growth rate k. Horizontal axis is time in 
units of inverse thermodynamical phase transition temperature, T c — 1 , and vertical axis is 
ln{[R(t) — R CT ]/[R(0) — -Rcr]} = y- Initial growth rate is given by k = lim t ^o ^jf (when 
R(0) -» R cr ). In the figure v = 1T C , £ cr = 1.07243T" 1 , R CI = 19.695TT 1 , and R(0)/R cr = 
1.01. 

bubble radius is 



where the dimensionality of space, d, is explicitly visible. In the case d = 3 
Eq. (H|) follows from this by comparing with Eq. (||) . Eq. (S) states clearly how 
in the real world the initial growth of bubbles is qualitatively different, much 
slower, than in one spatial dimension. 

The opposite approach is to let the bubble to expand in a hydrodynamical 
computer simulation, and to measure the initial growth rate numerically. Bub- 
ble radius is defined to be the distance where the tension or gradient energy 
has the maximum. The value of k can be read from Fig. ^ as the slope of the 
curve at origin. In the example case this gives k — (0.0052 ± 0.0001)T C , which 
coincides well with the analytical estimate from Eq. (||), n w 0.00516T C . 

These results can be applied in the analysis of thermal cosmological phase 
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transitions. In the case of relativistic heavy-ion collisions there is severe doubt 
on the validity of the general framework, nucleation in thermal systems. 
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